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5^ ■ Abstract 

In this paper we present a Priedmann-Robertson- Walker cosmological 
j_j ' model conformally coupled to a massive scalar field where the WKB approxi- 

mation fails to reproduce the exact solution to the Wheeler-DeWitt equation 
for large Universes. The breakdown of the WKB approximation follows the 
same pattern than in semiclassical physics of chaotic systems, and it is asso- 
ciated to the development of small scale structure in the wave function. This 
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result puts in doubt the "WKB interpretation" of Quantum Cosmology. 
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I. INTRODUCTION 

Quantum Cosmology |]l| tries to provide a complete description of the universe. Since the 
present universe obeys classical laws with great precision, an acceptable Quantum Cosmology 
must predict a "quantum to classical" transition at a certain time in the cosmic evolution. In 
the same manner, it must provide an interpretation of the initial conditions of the classical 
evolution, in terms of quantum processes. For these reasons, it is clear that the problem 
of the correspondence between classical and quantum behavior is a central point in the 
development of quantum cosmology. 

To explain our observations it is enough to interpret the wave function of the universe 
in the semiclassical limit [||. Moreover, it has been argued that it is sufficient to have a 
Wentzel-Kramers- Brillouin (WKB) state or a Gaussian state to have a valid semiclassical 
description 0. This is the basis of the so-called "WKB interpretation of Quantum Cos- 
mology" . According to this interpretation, we can understand the semiclassical limit as the 
limit rup — > oo, where rup stands for the Plank mass. In this limit the wave function of the 
universe, may be understood as a superposition of solutions of the WKB type 0], 

v|/ = ^C„exp[^5«] (1) 

where the Sa are solutions to the classical Hamilton Jacobi equation and the IC^p may be 
interpreted as the relative probability the a classical solution. However, we want to obtain 
the semiclassical limit when the universe has expanded beyond a certain size, and not in 
the limit mp -^ oo. The m^ -^ cxd limit of the wave function of the Universe would provide 
nevertheless information on the large size behavior, if the limits nip -^ oo and a -^ oo 
commuted, where a stands for the radius of the universe in a Robertson- Walker model. 
Although more recent papers have suggested the necessity of including decoherence to the 
WKB limit, in order to obtain a truly semiclassical state ||^, it is commonly assumed in 
the literature on quantum cosmology that the wave function of the universe adopts a WKB 
form at some point in the cosmic evolution, and retains it for ever after. 



This problem is similar to the problem of the commutativity of the limits ^ ^ and 
t -^ CX3 in ordinary quantum mechanics. The semiclassical wave function has been studied 
in great detail, and it has been shown that for irregular systems (chaotic systems) there 
exists a breakdown of the validity of the WKB approximation PJ^. The existence of chaotic 
cosmologies [§] suggests that a similar breakdown in the WKB approximation might occur 
in cosmological models. This would put in doubt the WKB interpretation of quantum 
cosmology. 

However, the problem of quantum cosmology cannot be referred directly to the quan- 
tum mechanical problem, because the Wheeler-DeWitt equation, the equation that governs 
quantum cosmology, is of second order. Therefore, we need to analyze the problem from 
first principles; it is convenient to begin such analysis with a simple, exactly solvable model. 

In this paper, we present such a model, namely a spatially flat Friedmann-Robertson- 
Walker (FRW) model coupled to a scalar field. Choosing a suitable coupling of the field to 
the radius a of the universe, we solve the Wheeler-DeWitt equation exactly. For the same 
model we calculate the semiclassical wave function, and we compare the evolution of this 
wave function with the exact one, checking the validity of the WKB approximation. 

The results obtained in this paper confirm the suspicion that the WKB approximation 
breaks down in quantum cosmology for large universes. This puts in doubt the WKB 
interpretation in its original formulation. The WKB approximation may still be valid, 
however, in more realistic models, for example, those including decoherence. 

II. THE MODEL 

Our cosmological model assumes a FRW spatially fiat geometry, conformally coupled to 
a real massive scalar field $. We shall suppose a ^(r, 77) = exp{i k.r)(f){ri)/a{ri), with A; << 1. 
This is a dynamical system with two degrees of freedom, a and 0, and Hamiltonian [§] 

H = \{^^ + vr,/,' + {k^ + mV(a)]02} = Q (2) 

Where •n^ y Pa are the momenta conjugated to and a respectively. 



For the numerical calculation we choose the potential 

oo 

V{a) = Y.{2n-l)e{\a\-5an)5a'^ (3) 

ra=l 

where 0(x) is the Heaviside step function and 5a is the step length (fig.|l|). 
Introducing dimensionless variables Pa/mpk"^ -^ Pa, anip/k — > a, T^,f,/k'^ -^ n^, (p/k -^ 0, 
m/nip -^ m, the adimensional Haniiltonian reads 

-| oo 

H = 7^{-Pa^ + ^/ + [1 + m^ ^(2n - 1) 0(|a| - San)5a^](l)^} = (4) 

^ n=l 

We can reduce the system's freedoms by one, using the Haniiltonian constraint. Hamil- 
ton's equations for the reduced system in the range n6a < a < {n + l)6a, are 

d(j) Sh n TT , . 

da dn -^/vr^ + Un^cj)'^ En 

da~ 6(j)~ 0f^T^J70^ ~ En ^ ^ 

Where now a plays the role of time, h = —Pa = En = y/'^'^ + ^n^'P'^ and 

uJn{a) = [l + m^nHa^Y'^ (7) 

is the frequency in each step. These equations of motion correspond to an harmonic 
oscillator of frequency ^, so the solution for the reduced system is 

0(")(5a) = r^' = cl>- cosi^6a) + ^ sin(^5a) (8) 

11/ UJn ^ 

vrW(5a) = Ti-^' = -,;„<^-sin(|^H +vr"cos(|^H (9) 

These equations define a stroboscopic map of period 6a, where 0" and Tr" are the values 
of the field and its conjugated momentum in the border of the nth step. 

As we cross the threshold from one step to the next, the energy changes from En to 
En+i- This non-conservation of the energy is a signal of the non-integr ability of this model. 



A. Semiclassical Wave Function 

Having found the solution to the classical problem we proceed to construct the semiclas- 
sical wave function as discussed by Berry |^. The idea is to associate wave functions \l/(0) 
to N-dimensional Lagrangian surfaces S, in the phase space ((f), tt^) (in our case, N=l). 
The association between E and \1/ is purely geometric, the dynamics of the problem will be 
introduced later, when we evolve the surface. We take the initial surface S^ as the invariant 
curve of the classical Hamiltonian for a massless field (m=0). 

h' = TT^ + 02 = E,' (10) 

Where Eo = 2u!o{k+l), Uo = 1, and k is a. natural number, as in Sommerfeld quantization 
rules. 

Our task is to associate S to a wave function 

^{<jy) = A{<jy)exp{zS{<j>)) (11) 

Where g 

A\^) = K\^\ = K\^\ = K\^\ (12) 

with K a constant of proportionality, and 5 is a solution of the classical Hamilton Jacobi 
equation, parametrized by /. 

If we now generalize this to a curve S„, obtained by evolving Sq n times through the 
classical map of equations (|^) and (^, the probability density associated to S„ is given by 

' t/0" ' ' d<p° " rf0" ' ° ' rf0" ' ^ ' 

To obtain the phase of the associated wave function we proceed as follows. The phase 
difference between two points separated 5(p on the surface E is 

ASo = So{(p°,I)-So{0,I)= I Ti°d(p°= f2Icos\^)d^ 



= I^ + ^ — - (14) 

^ 2sin(2(/?) ^ ^ 



And the phase difference between two points on the curve S„ (fig.||), reads 

A^„ = ASn-i + ^(vr/0" - 7r^"-V""' - EJa) (15) 

Where we neglected a global phase cr„(0"'(O), /). 

Since the function TT^{<f)) is multivalued, the wave function corresponding to the curve 
S„ is given by the superposition principle as 

vl/"(0) = E^""(0) exp[^5„"(0, J) + |/x] (16) 

where a labels the different branches of S„ for a given value of (p, and /x is the Maslov 



index associated to each fold of the curve S„ ||10| , p!l . 



B. Solution to the Wheeler-DeWitt Equation 

Following the canonical quantization procedure, we obtain the Wheeler-DeWitt equation 
for this model: 

^[^-^ + -.^(a)0^]^(a,0) = O (17) 

where uJn{a) is defined in eq.(^ and we have chosen the factor ordering so that the term 



in second derivatives becomes the Laplacian operator in the minisuperspace metric []T2 
Within the range n6a < a < {n + l)6a, we may expand 



oo 



^n{aA) = T.[^nFn' {a, <j)) + BJ FJ~ {a, <f))] (18) 



J 



where 



and 



Fj^ia, 0) = ^= exp[TVE„^a]$„^(0) (19) 



*n^'(0) = (v)'^'(2W^ exp[-ic.„0>^-(v/^0) (20) 



where hj{x) is the Hermite polynomial of grade j [|l^ and En'' = ujn{2j + 1 



Asking for continuity in the wave function and its normal derivative at the border of 
each step, we obtain the recurrence formulae for the coefficients of the expansion 

oo 

An+i' = Y.[Ar,^{Fn+^\FS)+Bj{Fn+f,Fn'')] (21) 

i 

oo 

5,+/ = ^[A„^(F„+/", F„^'^) + 5„^(F„+/', F„^")] (22) 

3 

Where {g^ f) is the Klein Gordon inner product 

(,.;) = ,:/,^(,.V_;^, (23) 

From the equations (|l^) to (^31) we finally obtain \l'n(a, 0) at the border of each step, 
from which we calculate |\l/n(a*, 0)p. This squared amplitude will be compared to that 
of the Semiclassical wave function. Also, we obtain the Klein Gordon charge, defined by 
the scalar product of eq.(p3|) as Qn = (\I^„(a, 0), \I^„(a, 0)). This charge is conserved by 
the Wheeler-DeWitt equation. We shall present the details of this recurrence relations in 
Appendix A. 

III. RESULTS 

Figures ^d and ^d show unfolding of the curve S evolved from an initial curve Sq through 
the classical map of eqs.(|D and (|^). The most outstanding feature is the development of 
spiral structures or "whorls", as Berry calls them 0. These are associated with invariant 



curves around a stable fixed point of h. They can arise, for example, in the twist map [jT4 
provided the angular frequency depends on the radius. In this case, points at different radii 
rotate around the central fixed point at different rates. Therefore radii map to spirals, and 
parts of S traveling close to stable fixed points will wrap around them, as it can be seen 
in figures ^ and ^. For these two figures (|^ and ^) we chose two different values of 5a, 
5a = 0.5 and 5a = 0.25 respectively. By the 10th step the "spiral galaxy" structure is clearly 
visible in the ^ curves. For the curves of figure |b, the spiral galaxy structure is visible 
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already at the 20th step (it must be noted that n = 10 and n = 20 correspond to the same 
"time" of evolution for the different figures), in spite of having a smoother potential. 

In figures ^a and ^ we can appreciate the semiclassical wave functions associated to the 
phase curves in figures ^ and ^d respectively. The most striking features are the caustic 
spikes, which proliferate as n increases and the classical curves curl over. 

The graphs of I^E'P, obtained from the exact solutions of the Wheeler-DeWitt equation, 
are shown in figures §3 and ^b for two different values of (5a, 5a = 0.5 and 6a = 0.25 
respectively. To study the corresponding quantum map (i.e. eq. |18D we must choose as 



initial state \E'o a wave function associated to Sq. Because Sq is an invariant curve of ho, 
\E'o must be an eigenstate of this operator. For the curves mapped on figures ^d and ||b, \E'o 
was chosen to be the 5th eigenstate of Hq. The graph of |\E'p for n = clearly shows the 
association with the initial curve Sq, with maxima at the caustics of the projection of S, and, 
between these, near harmonic oscillations. These oscillations correspond to the interference 
of two waves, associated with the intersection of 71^(0) with the fiber = constant, which 
make up the WKB function of figs, ^a and ^. 

By the 10th iteration (20th for 6a = 0.25) the WKB \^\^ (figs. |a (resp. |a)) has 
developed a complexity which cannot be observed in the real I^E'p. The caustics are evidently 
related to features of the wave function until n = 10 (resp. n = 20) but thereafter there is 
no obvious association. This is related to the fact that for large n the neighboring caustics in 
fig.Qa (resp. fig.^) are closer than the de Broglie wavelength, so they cannot be associated 
with features of \1' 0. This suggests, in the spirit of the smoothing procedure, that a better 
match between classical and quantal calculations will be obtained if we smooth, in some 
sense, the graph of I^E'p. Decoherence, or some other process that can eliminate details of 
the WKB function, could provide a suitable smoothing mechanism. 

In our numerical calculation of the exact wave function, we have computed only the first 
100 terms of the defining series eq. (p^ . As a check that this truncation does not impair the 
accuracy of our results, we show in figs. ^ and |^a a logarithmic plot of these coefficients. It 
is clearly seen from these plots that the coefficients decay exponentially, thus ensuring that 
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the tail of the series does not influence the value of the wave function within the accuracy of 
our calculations. Observe that the decay of the expansion coefficients for 5a = 0.25 (fig.^) 
is markedly faster than for 6a = 0.5 (fig.^a), as we should expect for a smoother evolution. 
As another check of the numerical accuracy, we computed the Klein Gordon charge for each 
step, verifying that it was conserved. 

IV. DISCUSSION 

In this paper we have presented a simple cosmological model where the WKB approxima- 
tion fails to reproduce the exact solution to the Wheeler-DeWitt equation for large Universes. 
The breakdown of the WKB approximation follows the same pattern than in semiclassical 
physics of chaotic systems, and it is associated to the development of small scale structure in 
the wave function. This result puts in doubt the so-called "WKB interpretacion of Quantum 
Cosmology", at least in its original formulation 0]. More complex models may provide a 
smoothing mechanism for the wave function, thus restoring the WKB approximation. 

In this model, the breakdown of the WKB approximation follows from the joint action 
of two effects, namely, the twist caused by the dynamics within each step, and the "particle 
creation" effect, that produces the excitation of higher eigenstates at each threshold. Since 
each of these effects are present in the continuum limit, 6a -^ 0, we may conclude that the 
WKB approximation will not be restored there. This is put into evidence by the fact that 
halving 6a leads to a stronger, rather than weaker, faillure of the WKB aproximation. 
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APPENDIX A: RECURRENCE FORMULAE FOR THE EXACT WAVE 

FUNCTION 



The scalar products in equations (|2^) and (|23|) may be computed from the recurrence 
relations 



'F _l/+ F ^' 

y-'- n+l ; -"■ n 



\^ n+l ) -i n ) 



(F„+/-,F„^'+) 



^ f:^'^J^li' exp[^(VK+i^ - ^EJ)a%n + l-t\n-j) (Al) 

Z{hn+1 i^n ) ' 



E., 



n+l 



-En 



2(E„+i*E„^)V4 



exp[z(VK+i* + \jEn')a*]{n + l-i\n-j) (A2) 



K+i* + \/^. 



2(E„+i*E„^)i/4 



exp[z(-VK+i* - VK^)«*](^ + i;^l^;j) (A3) 



{Fn 



+1 j-f^n 



^ n. 



En+1 



E„ 



2(F„+i^F„^)i/4 



f-—^exp[z{-^E^^^,^ + .jE^^)a*]{n + l;t\n;j) (A4) 



In these equations, 



{n + l;i\n; j) = J dx $n+i'*(x) ^J{x) 



These brackets can be written in recursive form 



(n + l]j\n;i) = -^\ — -. — {n + 1; j - 2\n; i) -\ \ -{n + 1; j - l\n] i - 1) 

Otn\ J Otn\ J 

for i < j, and 



/3„ /i — 1 1 / j 

{n + l;j\n;i) = \ — -. — {n + I; j \n; i - 2) -\ \ - {n + l;j - l\n;i - 1) 

an V « an V « 



if j < i, where |T3 



(n + l;0|n;0) 




and 



a^ 



1 r h^n+1 j ^n n 

2 Y u;„ Y u;„+i 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



Un is the step frequency given in eq.(j^) 



UJr. 



^n+l 



(AlO) 
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FIGURES 

FIG. 1. Comparison of the potential V{a) = c? with the potencial of the equation (^). 

FIG. 2. The diference of phase between two points on the curve S„ may be computed as the 
diference of phase on the curve Sn-i plus the phase gain from evolutin in ba 

FIG. 3. (a) Semiclassical Wave Function associated with the evolution of the curve S for a 
value of ba = 0.5. (b) Evolution of the curve S through the classical map, defined in equations @ 
and (P), for a value of ba = 0.5. 



FIG. 4. (a) Weight of the coefficients of the defining series of equation ( |T8|) for a value of 
ba = 0.5. (b) Evolution of the exact Wave Function of the Universe through the quantum map, 



defined in equation (18), for a value of (5a = 0.5. 



FIG. 5. (a) Semiclassical Wave Function associated with the evolution of the curve T, for a 
value of ba = 0.25. (b) Evolution of the curve S through the classical map, defined in equations 
(M and (P), for a value of 5a = 0.25. 



FIG. 6. (a) Weight of the coefficients of the defining series of equation ( [Tsl) for a value of 
ba = 0.25. (b) Evolution of the exact Wave Function of the Universe through the quantum map. 



defined in equation (18), for a value of (5a = 0.25. 
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